Influence of the magnetoelectric coupling on the electric field induced magnetization 
reversal in a composite non-strained multiferroic chain 
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We study theoretically the multiferroic dynamics in a composite one-dimensional system consisting 
of BaTiC>3 multiferroically coupled to an iron chain. The method treats the magnetization and the 
polarization as thermodynamic quantities describable via a combination of the Landau-Lifshits- 
Gilbert and the Ginzburg-Landau dynamics coupled via an additional term in the total free energy 
density. This term stems from the multiferroic interaction at the interface. For a wide range of 
strengths of this coupling we predict the possibility of obtaining a well-developed hysteresis in the 
ferromagnetic part of the system induced by an external electric field. The dependence of the reversal 
modes on the electric field frequency is also investigated and we predict a considerable stability of 
the magnetization reversal for frequencies in the range of 0.5 — 12 [GHz]. 

PACS numbers: 85.80.Jm, 75.78.-n, 77.80.Fm, 75.60.Ej, 77.80.Dj 
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I. INTRODUCTION 

Multiferroics (MFs), i.e. materials that possess ferro- 
magnetic and ferroelectric properties, have recently at- 
tracted significant research-^— as they hold the promise 
of qualitatively new device concepts such as electric field 
induced magnetization switching at low heat dissipa- 
tion. In addition to single-phase MFsi, i.e. materials 
with ferroelectric (FE) and ferromagnetic (FM) or an- 
tifcrromagnctic (AFM) ordering (e.g. BiFcOg^), com- 
posite MFs£ are in the focus of current research. These 
are systems realized as heterostructures of a wide range 
of different materials^. This diversity is hoped to com- 
pensates for the scarcity of single-phase MFs for room- 
temperature applications and to offer new routs to the 
engineering of the magnetoelectric coupling and its en- 
hancement, e.g. via an appropriate multilayer stacking. 
A typical example for such composite FE/FM junctions 
is BaTiOa/Fe that was predicted^ and studied^ theoret- 
ically. Experimentally BaTiOs/Fe was successfully real- 
ized and characterized ! 10 i 11 In addition, a recent publi- 
cation reported on a considerable coupling between the 
iron and the barium titanate which renders possible a 
change of the magnetization of the iron layer when an 
electric field is applied, even at a room temperature^ 
Another important issue concerns the significant varia- 
tion of the tunnel magnetoresistance depending on the 
polarization direction of BaTiC>3 layerii^ All these ex- 
perimental findings are very promising steps towards the 
creation of memory elements^ 3 - based on composite MFs, 
which could be written by electric field pulses. Prob- 
lems appearing on the way to achieve this goal include 
the optimization of the geometry of FE/FM layers and 
the improvement of the interaction between them. This 
paper is dedicated to the latter issue. 

In contrast to our earlier studies^ - — where we focused 
on a more simple tetragonal phase of barium titanate in 
which the pcrovskite exists at room temperature^ here 



we consider the rhombohedral phased at zero Kelvin 
with the corresponding ferroelectric potential. Under 
these circumstances it is possible to neglect the thermal 
stochastic effects on the effective fields. To improve the 
accuracy of the model, we also included the long-range 
interactions for both the FE and the FM materials and 
the effect of the FE depolarizing field. 



II. THEORETICAL FORMULATION 

We consider a chain composed of FE (BaTiOs) and 
FM (Fe) sites as shown in Fig. 1. The aim is to describe 
the evolution of the magnetization (polarization) under 
the influence of external fields. We adopt a Ginzburg- 
Landau phenomenology based on coarse grained order 
parameters that formally result from an averaging of 
the relevant microscopic quantities over an appropri- 
ate cell. These cells are in our calculations for the 
whole chain cubes (called hereafter sites) of equal vol- 
ume Op E = cip M = a 3 = 5x5x5 [nm 3 ]. We will show 
results for a chain formed by five FE and five FM sites. 
Calculations performed for longer chains^ showed that 
the total magnetization reversal is hardly achievable for 
a chain that has more than ten magnetic sites. 

The coarse grained total free energ y 14 : 15 reads 

Ftot = a 3 FpE + a 3 FpM + Econ- (1) 
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The depolarizing energy density Pdep— - (cf. Fig. 1) 
reads 
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which involves the dead layer permittivity £dl, the di- 
electric constant in vacuum eo, the metal screening length 
Am, the FE cell size a, and the cell number N. The ferro- 
electric nearest neighbors coupling energy density F^p L 
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where kfe is the ferroelectric coupling constant. The 
dipole-dipole interactions Ppf^ is 
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Here £fe is the ferroelectric permittivity, is a unit 
vector parallel to the line joining the centers of the dipoles 
Pi and Pfe and Ui k is the distance (measured in units of 
a) between the two dipoles. 

The energy term Pext stems from the applied electric 
field E 



Pext = -^P-P- 



(6) 



Analogously the ferromagnetic part is characterized by 
the coarse grained free energy density 
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which consists of the (uniaxial) magnetocrystallinc 
anisotropy contribution^ 
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with the anisotropy constant Ki and the saturation mag- 
netization Ms- The nearest- neighbor exchange interac- 
tion has the form22 
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where A is the interaction constant, and the dipole-dipole 
interaction is 
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where /io is the susceptibility constant. The coupling 
between the ferroelectric and the ferromagnetic parts^ 3 - is 
due to the mobile spin-polarized electrons accumulated at 
the interface in order to screen the electric polarization in 
the FE part^ A change of the accumulated spin density 
(e.g. due to a change of the electric polarization) will act 
with a torque on the magnetization. This is however a 
surface effect restricted to the region in the vicinity of 
the interface (which is the reason why we arc considering 
short chains). In other words, only surface cells (those 
with index 1 in Fig. 1) will participate in this coupling 2 ^ 
which contributes with the interaction energy 
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The time dynamics of the polarization Pj and the mag- 
netization Mj of the individual sites are obtained by 
propagating the Landau-Khalatnikov (LKh) 26 ' 27 equa- 
tion 
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and the Landau-Lifshitz-Gilbert (LLG) ' equation 
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The coefficients entering into the LKh are the viscosity 

constant 7„ and the effective electric field Pfej- The LLG 

equation involves the gyromagnetic ratio 7, the Gilbert 

damping coefficient otm, the saturation magnetization 

-* 1 <5Ptot 

Ms, and the effective field Pfm> = 7 — ■ 

a 3 6Mj 

From a computational point of view the appropriate 
choice of the cell size a is important. If it is cho- 
sen too small, the coarse graining procedure to obtain 
the macroscopic quantities Pi and Mj becomes ques- 
tionable and one faces in addition problems with the 
supcrpara- magnetic/electric limits. For a cell too large, 
a multi-domain state sets in. The cell size used in our 
calculation is therefore a = 5 [nm] for both FE and 
FM parts of the chain. The material parameters were 
set as those of barium titanate at T = [K]. Specif- 
ically, we choose a FE i = -1.275 x 10 s [V-m/C]^ 
f3 FE1 = -2.045 x 10 9 [V-m 5 /C 3 ]^ /3 FE2 = 3.230 x 
10 s [V-m 5 /C 3 ]^ 7FE1 = 9.384 x 10 9 [V-m^C 5 ]^ 
7FE2 = 4.470 x 1 " 9 rv--.9/r<5i.30 
10 9 [V-n^/C 5 ]^ 



4.470 x 10 9 [V-m 9 /C 5 ]^ 7fe 3 = 4.919 x 



7„ = 1.5 x 10 



" 5 [V-m-s/C]^- P s 



0.499 [C/m 2 ]^K FE = 2.04 x 10 s [V-m/C]^. e FE = 
1643 -Fjjg depolarizing energy density Pdep includes 
the permittivity of the metallic electrodes, or the so 
called "dead layer" j2i which can be lower than the per- 
mittivity of the material. This parameter is difficult to 
measure 3 ^. In the numerical calculation we used the rea- 
sonable value £dl = £fe/2 = 82. The material of the 
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FIG. 1. Schematics of the composite multiferroic structure 
formed of five ferroelectric and five ferromagnetic sites. The 
initial state for both the FE polarization and the magnetiza- 
tion is chosen as random. 



FM layer is taken as iron with the following parameters 
at T = [K]: a FM = 0.5^ K x = 4.8 x 10 4 [J/m]^2 
M s = 1.71 x 10 6 [A/m]^2 A = 2.1 x 1CT 11 [J/m}M 

A special role in the multiferroic dynamics is played 
by the coupling strength A. We estimate the strength 
of the magnetoelectric coupling parameter related to 
the BaTiOa/Fe- interface using the proposed ab-initio 
expressions^ as = f-oAM/Ec, where the surface mag- 
netoelectric coupling as is defined as the change of the 
surface magnetization AM for the electric coercive field 
Eq. The coercive field can be estimated as Eq w 
-Ps/(£fe£o)- Keeping in mind that the coupling energy 
(|11[) might also be expressed through the induced mag- 
netization AM$ at the interface and the net ferromag- 
netic magnetization M as Eqon = — J/(M|)AMi • M, 
we finally obtain A = Jas/i^oM^a £fe£o)- For as = 
0.2 • 10" 17 [T m 2 /V}2 the last expression yields A « 
0.063 [s/F]. This numerical value did not result in any siz- 
able interaction between FE and FM layers in the model 
considered. We attribute this situation to the difficulty 
in the numerical definition of A, because the expression 
for the coupling constant depends strongly on the size of 
the interface cell a. Under these circumstances, the most 
adequate approach was to consider A as a variable, with 
the aim to obtain the general picture of magnetization 
reversal induced by an electric field as a function of the 
magnetoelectric coupling. 

As in our previous studie a 14 ' 15 the MF-chain will be 
characterized by the averaged total polarization p A = 
{NPs)~ 1 ^2 i Pi and the averaged net magnetization 



To reverse the FM part via the FE part of the chain, 
we apply the harmonic electric field E x (t) = E sm(uit) 
with the amplitude Eo = 8x 10 7 [V/m]. As a frequency 
we choose at first lu/(2tt) = 2 [GHz]. The time profiles 
of the field reversal for the ferroelectric and the ferro- 
magnetic parts of the system are given in Fig. 2. As 
one can see, the ferroelectric part re-polarizes quickly 
within approximately 300 [fs] . The reversal modifies only 
the x-component of the polarization, with a very minor 
variation of the z-component in the vicinity of the point 
where p A changes its sign. This happens because under 
a negative bias the polarization vector of the pre-contact 
diminishes in magnitude to zero, and then reverses the 
direction along the field. The remaining cells flip their 
polarization vectors along the x-axis without a preces- 
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FIG. 2. The difference in the time scale for the ferroelectric 
(upper panel) and the ferromagnetic (lower panel) reversals. 
The FE/FM coupling constant is taken as A = 10 [s/F]. 



sion, causing a fast and a well-developed reversal. On 
the contrary, the magnetization reversal requires a much 
larger time up to a nanosecond. As one can see from 
the plot, the re-orientation of the magnetization vector 
involves a heavy precession with considerable deviations 
of the rriy and components from zero. These oscil- 
lations have a higher frequency at the beginning of the 
reversal process, gradually lowering the frequency and 
the amplitude till the equilibrium state is reached. Natu- 
rally, such a distinct time scale of the system components 
poses a significant problem for a proper modeling. To 
ensure the accurate numerical solution, one should keep 
the integration time step at femtoseconds, which drasti- 
cally increases the number of integration steps needed to 
follow the system dynamics on time scales adequate to 
observe the hysteresis curves formed under the field vari- 
ation with GHz frequencies. This large number of steps 
definitely will be an issue for the calculations of large sys- 
tems composed of hundreds of particles due to the need 
to evaluate the long range dipolc-dipole interaction fields 
in the ferroelectric and the ferromagnetic parts. In our 
case, the LKh and the LLG equations were integrated 
with the Heun method using time step At = 2 [fs] . 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 



To study the influence of the coupling constant on the 
behavior of the system, it is necessary to select param- 
eters that allow an easy and a reliable characterization 
of the hysteresis curve. We propose to use the values of 
the averaged polarization and the magnetization (p A and 
m A , respectively) at a zero applied field, as well as the 



4 



(a) 

























Hysteresis curve components: 

1 



■ p A , m A — p A , m A — p A , m A 

1 x 7 x y y z z 







-1- 



(b) 





























-1 



(c) 

























(d) 

























:(e) 












....Mm 


■ 1 1 1 1 1 ii i 



-100-50 50 100 -100-50 50 100 -100-50 50 100 -100-50 50 100 -100-50 50 100 



-1- 




-1 - 









k 


■\ 

/ 














i 


j 




1 1 1 1 1 1 1 11 


■ ■ ■ . i . . . . 



-1- 



-1- 



1 






— V 



-100-50 50 100 -100-50 50 100 -100-50 50 100 -100-50 50 100 -100-50 50 100 
E[MV/m] E [MV/m] E [MV/m] E [MV/m] E [MV/m] 

1.2 





0.6- 






;d 


- 


<oi 






-0.6 - 




-1.2- 




10 




1.2 -| 




0.6- 








- 






B 


-0.6 - 




-1 ? - 






100 -n 


a 


50^ 


> 






0^ 




-50: 


w 





'- • •■ 



10"- 



10" 



10° 
I [s/F] 



-100 

1 

100 
^ 50 

a o 

§o -50 

w 

-100 



10' 



10 2 




A, [s/F] 



FIG. 3. The dependence of the multiferroic reversal on the strength of FE/FM interaction A. The hysteresis curves are 
presented for: a) A = 0.1 [s/F], b) A = 0.255 [s/F], c) A = 2 [s/F], d) A = 12 [s/F], e) A = 45.5 [s/F]. The hystograms in 
the lower part of the figure depict the remanence (p^, mj^) end coercitivity (E^ E , Eq M ) for the averaged hysteresis curves 
obtained for ferroelectric and ferromagnetic parts of the chain. The applied electric field is characterized by the amplitude Eo 
= 80 [MV/m] and frequency oj/(2tt) = 2 [GHz]. 



field values Eq for which the magnetization or the polar- 
ization are equal to zero. If the hysteresis loop is properly 
formed, these points will define the remanence and the 
coercitivity of the hysteresis loop. If the hysteresis is ab- 
sent, it will become immediately notable, for example, by 
the existence of a single remanence point or the presence 
of various zeros of magnetization / polarization curve. 
To collect the data shown in Fig. 3 we integrated the 
system for 15 full field cycles and plotted the density of 
p A = p A {E = 0), m A = m A (E = 0) and £g E , E™ 
as a function of the coupling parameter A. These plots 
are shown in the bottom part of Fig. 3. This approach 
allows to sec if the hysteresis curve is well-defined and 
repeatable, which will result in sharp and dark lines on 
the corresponding plot. If a hysteresis loop is unstable 
and varies from cycle to cycle, one obtains the grayish 



bands instead of the thin dark curves. The light-colored 
noise for the positive fields on the coercitivity diagrams 
correspond to the points when the polarization / magne- 
tization changes sign before reaching the stable hysteresis 
curve. In some cases, such "self- adjustment" of the phase 
trajectory takes place during the several first cycles of 
the electric field, which complicates the definition of an 
exact threshold that would allow to remove this noise. 
On the other hand, as we are studying steady hysteresis 
curves observed during the 15 th field cycle, such initial 
noise is irrelevant and can be easily neglected. The most 
characteristic cases with FE and FM hysteresis curves il- 
lustrated in the upper part of the figure are marked with 
arrows in remanence and coercitivity diagrams. 

As seen from the figure the coupling strength A < 0.8 
[s/F] is insufficient to reverse the magnetization. For the 



5 



small value A = 0.1 [s/F] (Fig. 3a) the value of is con- 
stant and matches the saturation magnetization of the 
material, while the ferroelectric part manages to reach 
an entire hysteresis cycle. The steps at the edges of the 
hysteresis are caused by the discussed peculiarity of FE 
reversal mechanism that involves a vanishing of the first 
site before the field achieves the reversal value. For a 
stronger coupling A = 0.255 [s/F] (Fig. 3b) the torque 
rendered by the polarization "kicks" the magnetization 
with a frequency that is similar to the FM precession 
frequency, resulting in a periodic variation of m^, 
and m^. These oscillation modes may be interesting as 
a way to achieve a resonant magnetization precession. 
However, they do not allow to reach a saturation magne- 
tization and to form the hysteresis loop, which is required 
for the operation as memory devices. 

The full ferromagnetic hysteresis loop emerges for A > 
0.8 [s/F], quickly reaching the full saturation value Ms- 
Due to the fact that the FM part is much slower than 
FE, the magnetic hysteresis has a larger coercitivity, how- 
ever. This is most clearly seen for the case of a relatively 
weak coupling A = 2 [s/F] (Fig. 3c), where FM dynam- 
ics definitely can not follow FE dynamics at the desired 
speed. We observe pronounced oscillations in the my 
and components upon the field-induced reversal. It 
is worth mentioning that the FE part of the system is 
not influenced much by the variation of the magnetiza- 
tion for these values of A, for the ferroelectric hystere- 
sis exhibits only a minor enhancement of the remanence 
and has almost the same coercitivity. An increase in the 
magneto-coupling strenght A to 12 [s/F] (Fig. 3d) seems 
to yield desirable results: the magnetization hysteresis 
is saturated and narrow, with a fast relaxation of the 
rriy and components after the full reversal of m*. 
As the FE/FM feedback becomes more pronounced, the 
hysteresis of FE shows a larger coercitivity and rema- 
nence. The observed "mirror-symmetry" of the FE and 
FM hysteresis loops has its origin in the opposite sign of 
the coupling terms in the LKh and LLG equations. For 
a large coupling A = 45.5 [s/F] (Fig. 3e) both hystere- 
sis curves deteriorate, as the ferromagnetic part starts to 
"hold" the FE sites, hindering their reversal until a con- 
siderable electric field is applied. A further increase of 
the coupling strength destroys the hysteresis of both the 
FE and FM layers as A exceeds 60 [s/F]. 

To address the frequency dependence of the composite 
multiferroic reversal, we calculated the remanence and 
coercitivity diagrams (Fig. 4) for the frequency range 
ui/(2ir) =0.5-12 [GHz]. As one can see from the figure, 
for a low frequency us/(2ir) < 7.2 [GHz] the system is 
well-tuned, featuring a complete magnetization reversal 
triggered by the ferroelectric component of the struc- 
ture. The increase in the frequency results in a linear 
increase of the coercetivity field for the ferromagnetic 
hysteresis. The low-frequency hysteresis has a sharp and 
a well-defined shape (Fig. 4a). The slow response of 
the FM to the faster oscillations of the FE polarization 
becomes noticeable under the increase of oj, which man- 



ifests itself by rounded "corners" of the hysteresis curve 
(Fig. 4b) that becomes more pronounced for a frequency 
around 7 [GHz]. At this value, one can also observe a 
slight decrease in the magnetic remanence (Fig. 4c). 
When to/ (2ir) exceeds the characteristic frequency of iron 
27A'i/[(l+a| M )M s /i ] = 7.89 [GHz], the behavior of the 
system changes drastically. The FM part does not follow 
the fast FE dynamics so that the curve of the magnetic 
hysteresis features several "breaks" and its coercitivity 
slightly drops down. It is essential to stress that, because 
we are dealing with a complex multi-site system driven 
by FE/FM interaction, the changes of the magnetiza- 
tion dynamics are noticeable for the frequency of electric 
fields slightly lower than the characteristic frequency of 
iron. The magnetization reversal is accompanied with a 
strong deviation of m y and m z components (Fig. 4d). 
A further increase of the frequency destabilizes the sys- 
tem even more, with the magnetization hysteresis looking 
"skewed" and the coercitivity values "jumping" back and 
forth between two branches (Fig. 4e) . It should be men- 
tioned that the FM remanence lowers constantly for a fre- 
quency u>/(2ir) > 6.5 [GHz], signaling that the achieved 
operation modes are not very promising for applications. 
Also it is important to stress that the considered value 
of the coupling constant A = 12 [s/F], while allowing a 
good influence of the FE over the FM layer, provides a 
minimal feedback that can be seen in a lack of the dras- 
tic changes of the remanence and the coercitivity for the 
FE part of the system. It goes without saying that such 
a pronounced unidirectional connection is an important 
feature for a device that is aimed to control the magnetic 
dynamics with an electric field via FE/FM coupling. 
Based on these simulations, we would recommend to use 
the described compound multiferroic system under low 
GHz frequency of the applied electric fields. 



IV. CONCLUSIONS 

We reported on a full-scale magnetization reversal in a 
composite multiferroic chain using an applied harmonic 
electric field with frequencies 0.5 — 12 [GHz]. The dynam- 
ics and in particular the reversal depend sensitively on 
the strength of the magneto-electric coupling A. Hence, 
tracing the dynamics should deliver information on this 
coupling. For a weak coupling A < 0.8 [s/F], the ferro- 
electric part shows well-defined hysteresis loop which may 
result in periodic oscillations of the ferromagnetic part. 
A saturation magnetization Ms may not be reached how- 
ever. A magnetic hysteresis opens when A grows above 
0.8 [s/F]. For a coupling constant strength in the ranges 
10 — 20 [s/F] , the multiferroic system has an optimal per- 
formance with both FE and FM hysteresis curves featur- 
ing a high remanence and a low coercitivity. It is im- 
portant to highlight that the range of A that corresponds 
to the FM hysteresis with the most definite shape does 
not vary with the frequency of the applied electric field. 
When the coupling constant becomes too strong (exceed- 
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FIG. 4. The frequency dependence of the FE/FM reversal for A = 12 [s/F]. The most characteristic hysteresis curves are given 
for the following values of the field frequency w/(2tt): a) 0.5 [GHz], b) 4 [GHz], c) 7 [GHz], d) 9 [GHz] and e) 11.4 [GHz]. 



ing 40 [s/F]) the system does not show any hysteresis. 
This degradation is caused by the fact that the ferroelec- 
tric dynamics is strongly disturbed by the ferromagnetic 
layer. 

We considered here one type of magnetoelectric cou- 
pling based on a screening model at the interface of the 
FE/FM parts. A different type of magnetoelectric cou- 
pling mechanism, e.g. a stress-strain coupling or a dy- 
namical Dzyaloshinskii-Moriya coupling^, may qualita- 
tively alter the coupled dynamics (even when the cou- 
pling strengths are comparable), for the corresponding 
functional forms entering the equation of motion are dif- 



ferent in general. Hence, in addition to the relevance for 
application, an attractive feature of studying the com- 
posite dynamics is that it may deliver some details on 
the underlying multiferroic coupling mechanisms. Stud- 
ies along this line are currently underway. 
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